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I  his  paper  describes  an  entrainment  model  for  cohesive  sediments  that  is  based  on 
a  power-law  expression  for  the  excess  shear  stress  and  a  total  entrainment  coefficient. 
Total  entrainment  includes  terms  for  consolidation,  bioturbaiion,  and  all  other  processes 
i base  entrainment  coefficient).  The  model  is  used  to  predict  entrainment  rates  for 
cohesive  sediments  from  Lake  Erie,  the  Tamar  Estuary,  Long  Island  Sound,  and  the  Fox 
River,  Wisconsin.  The  base  entrainment  coefficient,  which  is  estimated  using  samples 
with  the  least  post-deposit  ional  modification,  is  unique  for  each  sediment  suite  because 
a  includes  environmentally  sensitive  processes  like  mineralogy,  salinity,  organic  carbon 
content,  etc.  Based  on  available  entrainment  measurements,  expressions  are  presented 
fur  the  con  soli  da  lion  and  bioturbation  coefficients.  The  model  is  evaluated  with  en¬ 
trainment  data  for  identical  sediment  that  has  been  either  consolidated  or  bioturbated 
and  the  comparison  is  encouraging.  A  comparison  of  predicted  and  measured  entrain¬ 
ment  rates  for  undisturbed  sediment  is  less  favorable  because  of  its  unknown  post* 
depositions!  history. 


1.  INTRODUCTION 

Cohesive  sediments  are  of  interest  because  of  their  role  in  shoaling  processes  in  cs~ 
tuaries  and  harbors,  contaminant  transport,  and  marine  structural  design.  The  primary 
mot  i  vat  ton  of  the  US,  Navy,  however,  is  their  impact  on  the  physical  and  acoustic 
properties  of  the  seabed  and  on  the  optical  properties  of  the  water  column.  The  physical 
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and  chemical  environment  determines  the  dynamics  of  cohesive  sediment  but  interac¬ 
tions  between  processes  within  the  seafloor  and  the  water  column  are  poorly  known. 
This  problem  is  solved  through  a  dependence  on  field  and  laboratory  measurements  m 
non-military  applications*  but  this  approach  is  generally  unavailable  for  military  opera¬ 
tions.  Thus,  more  general  methods  need  to  be  developed.  This  paper  discusses  one  such 
approach  that  relies  less  on  direct  observations  of  cohesive  sediment  behavior. 


2.  BACKGROUND 

A  substantial  research  effort  has  gone  into  understanding  erosion  (also  called  en¬ 
trainment  or  resuspension)*  deposition,  and  consolidation  Parlheniades,  1993: 

Tuorman,  1996)  of  cohesive  sediments  The  entrainment  of  fine  grained  sediment  is 
usually  parameterized  using  a  power- law  function  that  incorporates  the  bottom  shear 
stress.  For  example*  Ariathurai  ei  at.  (1983)  proposed  the  following  expression  for  the 
entrainment  rate*  E: 


CD 


where:  M  =  the  rate  constant  (kg/m2/s)  at  xh  =  2t*;  Tb  =  bottom  stress;  and  rc  =  critical 
shear  stress  for  entrainment*  which  has  been  found  to  increase  with  depth  for 
undisturbed  sediment  (a*.,  Sanford  and  Maa,  2001).  Note  that  in  this  paper,  we  will  use 
a  constant  value  of  0.04  Pa  for  %  when  it  is  not  available.  Partheniades  (1965) 
suggested  that  erosion  rates  for  mud  (60%  day),  while  being  constant  for  some 
experimental  conditions*  could  only  be  described  by  a  complex  function  with  a  large 
number  of  coefficients.  Parchure  and  Mehta  (1985)  found  that  the  erosion  rate  of 
deposited  beds  was  not  constant  in  part  because  of  consolidation  and  dewatering.  They 
turned  lo  the  theory  of  chemical  rate  processes  to  explain  the  dynamics  of  cohesive 
sediments,  thereby  implicitly  making  entrainment  dependent  on  inter-particle  bonds  and 
thus  the  bio-gcochemical  environment.  The  parameters  in  their  proposed  formulation 
were  determined  by  fitting  experimental  data.  Lavelle  el  al  (1984)  examined  the  data 
from  a  number  of  studies  and  proposed  the  following  expression  lor  the  erosion  rate. 


E  =  a]Tb!n 


(kg/m  Vs) 


a> 
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where  the  parameters  a  and  T|  were  determined  by  least -squares  fit  to  the  data  from  each 
set  of  experiments.  The  parameter  a  varies  by  three  orders  of  magnitude.  More  recently. 
Lick  et  at  (1995)  proposed  the  following  expression  for  the  resuspension  potential  for 
cohesive  sediments; 


(3) 


where  n,  n,  and  m  are  constants,  and  tj  is  the  time  after  deposition  in  days.  If  the  shear 
stress  is  held  constant,  the  entrainment  rate  can  be  found  from  sfTr  where  T  is  the  time 
interval  over  which  the  entrainment  remains  constant  (equilibrium  time  step). 
Ex  peri  mental  work  has  indicated  that  the  mean  concentration  of  entrained  sediment  in 
laboratory  flumes  reaches  a  steady  state  m  approximately  l  hr  (Lick,  1982).  The 
entrainment  rate  can  be  deduced  from  the  initial  slope  of  the  concentration  lime  series, 
implying  that  E  will  go  to  zero  when  an  equilibrium  concentration  is  reached.  Thus, 
some  data  arc  reported  as  a  mean  entrainment  rate  (or  first- hour  rate)  for  this 
equilibrium  interval  (see  Lick,  1982).  For  these  data,  wc  calculate  E  from  €  using  T  = 
3600  s.  Note  that  only  Eq,  1  is  dimensionally  correct.  This  brief  summary  is  intended  to 
show  that  there  is  no  unifying  theory  of  the  physical  processes  that  control  the 
entrainment  of  cohesive  sediments  and  that  most  entrainment  models  use  empirical 
parameters  that  are  derived  from  observations  of  sediments  from  the  area  of  interest. 

The  influence  of  biological  and  physical  processes  on  the  erosion  of  mudflats  has 
been  examined  in  the  INTRMUD  project  (Black  et  al.y  1998),  which  demonstrated  the 
importance  of  biostabilization  on  entrainment  (de  Brouwer  et  at ,,  2000;  Droppo  et  al, 
2001),  The  LISP  project  (Littoral  Investigation  of  Sediment  Properties)  (Daborn,  1991) 
also  examined  the  complex  interaction  of  microflora  and  fauna  in  the  intertidal  envi¬ 
ronment.  The  biostabilization  effect  of  bio  films  is  opposed  by  the  destabilization  caused 
by  blot  ur  bad  on  (Grant  and  Daborn,  1994;  Green  et  a/.,  2002).  Laboratory  work  with 
natural  sediments  has  shown  that  cohesive  sediment  credibility  increases  rapidly  with 
the  activity  of  infauna,  presumably  because  of  their  destruction  ol  inter* particle  bonds 
and  primary  depositional  fabric  (Tsai  and  Lick,  1988;  Limern  et  ai,  2002),  The  com¬ 
plex  interaction  of  these  biological  processes  w-ith  physical  processes  like  hydrodynam¬ 
ics  and  consolidation  necessitates  the  use  of  more  realistic  entrainment  models,  and 
carefully  planned  measurement  programs  (Tothurst  et  al .,  2000;  Lick  et  al,  1998). 
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3.  METHODS 

The  purpose  of  this  study  is  to  evaluate  sediment  entrainment  using  a  modular  ap¬ 
proach,  which  isolates  the  inherent  physical  propen ies  of  sediment  from  consolidation 
and  biomrbatfon  effects*  The  proposed  model  is  first  presented  and  the  method  of 
determining  expressions  for  the  coefficients  for  bioturbatkm  A#  and  consolidation  Ac  is 
then  discussed. 

3*1*  The  entrainment  model 

The  entrainment  expressions  discussed  in  the  literature  are  predominantly  derived 
by  fining  a  power-law  function  to  a  series  of  data  points  on  a  plot  of  xb  versus  £,  The 
power  varies  by  many  orders  of  magnitude*  depending  on  the  shear  stress  function.  For 
example,  powers  as  low  as  10 4  were  used  by  Lovelle  a  at  (1984)  in  Eq.  2,  whereas  a 
power  of  1  is  implied  by  a  linear  model  (Eq.  1)*  and  the  resuspension  potential  (Eq,  3) 
uses  a  coefficient  greater  than  L  These  expressions  all  have  leading  coefficients  that  are 
either  constant  or  variable.  This  simple  function  fits  a  range  of  data  because  the  sedi¬ 
ment  concentrations  of  interest  to  most  problems  are  small,  which  is  the  appropriate 
application  of  a  power-law  function.  This  function  ts  inappropriate  for  fluid  mud, 
however,  and  more  complex  models  are  required  (c.#.,  Krancnburg  and  Winterwerp, 
1997)* 

I  he  model  we  propose  to  describe  cohesive  sediment  entrainment  is  also  based  on 
a  power-law  expression,  but  physical  and  biological  processes  are  isolated  so  that  they 
can  be  evaluated  separately.  This  is  done  conceptually  by  decomposing  the  leading 
constant  into  independent  coefficients: 


where  Ac,  and  are  parameters  to  be  determined  from  observations.  This 
formulation  is  similar  to  Eq*  l  but  the  excess  shear  stress  power  m  can  vary  from  unity. 
The  parameter  Aq  is  analogous  to  the  leading  coefficients  used  in  other  expressions.  It  is 
treated  as  a  fundamental  entrainment  rate  for  specific  sediment,  which  has  undergone  no 
post-deposit  ion  al  modification  (r*c,*  no  consolidation  or  bioturbation).  The  parameters 
and  Ar  represent  biological  reworking  and  consolidation  processes,  respectively:  if 
no  consolidation  or  bioturbation  has  occurred,  both  AB  and  At  will  equal  l,  indicating 
no  effects.  We  neglect  biombihzation  biofilms)  Thus,  our  model  incorporates  a 
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fundamental  entrainment  rate  that  is  a  function  of  excess  shear  sirens,  and  which  is 
modified  by  post-deposit  ion  a  l  consolidation  and  bioturbation.  1  he  values  ot  the 
parameters  will  be  discussed  in  the  results  section. 

This  study  focuses  on  sediment  entrainment  and,  therefore*  suspended  sediment 
profiles  will  not  be  calculated  or  discussed.  Equation  4  will  be  evaluated  using  entrain¬ 
ment  rates  derived  from  annular  flumes  (Fukuda  and  Lick,  1980;  Tsai  and  Lick,  1988), 
the  Instrument  for  measuring  Shear  stress  In  Situ  (ISIS;  Limern  &  a/.,  2002)  and  di¬ 
rectly  measured  by  SEDFLUME  (McNeil  er  al ,  1996).  The  shear  stress  Xh  can  be 
calculated  from  the  hydraulic  flow  in  all  of  these  devices  and  thus  the  critical  shear 
stress  tc  can  be  estimated.  However,  the  definition  of  tc  is  partially  subjective  and  not  all 
workers  define  if  the  same.  For  example,  McNeil  et  ai  (1996)  define  Tc  as  the  shear 
stress  that  produces  an  erosion  rate  (measured  by  SEDFLUME)  between  10  and  10  * 
m/s;  thus  the  mass  erosion  rate  (kg/m7s)  is  dependent  on  the  sediment  bulk  density.  We 
use  measured  bulk  densities  to  calculate  entrainment  rates  for  comparison  to  the  model 
and  other  results.  Lmiem  et  nL  (2002),  however,  define  ^  as  the  stress  at  which  the 
entrainment  rate  increases  from  a  previous  steady  value  (see  their  Figure  3),  The  defini¬ 
tion  of  the  critical  shear  stress  reported  by  different  authors  is  dependent  on  measure¬ 
ment  methods,  and  this  variability  will  contribute  to  uncertainty  in  the  results  of  this 
study. 

3,2.  Experimental  procedure 

The  coefficients  A0.  Ac  and  Ae  in  Eq.  4  must  be  estimated  in  order  to  isolate  the  im- 
pacts  of  bioturbation  and  consolidation  from  other  processes.  Ideally,  this  would  be 
accomplished  using  entrainment  data  for  undisturbed  sediments,  but  ihis  approach 
cannot  be  used  because  the  post-depositional  history  of  such  sediments  is  unknown.  We 
will  use  homogeneous  sediment  that  was  deposited  from  slurries  to  find  functions  of  Ac 
and  Aa,  whereas  A0  will  be  independently  determined  for  each  sediment  suite.  The  data 
consist  of  measured  entrainment  rates  £,  the  resuspension  potential  £,  and  the  equilib¬ 
rium  concentration  Cm  at  different  shear  stresses,  water  contents,  and  consolidation 
intervals.  The  sediments  vary  somewhat  and  caution  is  necessary' m  developing  parame- 
ten  zat  ions  for  general  applicability.  Nevertheless,  this  first -order  calculation  of  the 
coefficients  will  assist  in  future  efforts  to  understand  these  processes.  For  ease  ol 
comparison*  all  entrainment  measurements  wall  be  transformed  into  an  entrainment  rate 
&  (kg/rrr/s),  which  requires  approximation  of  umeported  experimental  parameters  for 
some  data. 

The  base  entrainment  parameter  Ao  is  strongly  dependent  on  clay  mineralogy,  salin¬ 
ity,  and  seabed  chemistry.  These  variables  are  not  considered  in  this  study  and  thus  they 
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i hey  are  implicitly  included  in  A#,  For  this  reason.  At,  is  expected  to  have  a  unique  value 
lor  each  sediment  suite  discussed  in  this  report.  This  coefficient  should  be  determined 
Imm  sediment  ihal  has  undergone  minimal  posl-deposii tonal  modification.  For  each 
suiie  of  daia  {e.g. .  Lake  Erie  sediment  from  I  able  I ).  samples  that  are  assumed  to  be  the 
least  modified  by  posi-dcpositional  processes  are  used  to  fit  Eq.  4  and  thus  find  a  value 
of  40.  This  approach  assumes  that  Aa  =  Ac  =  I,  which  is  reasonable  for  newly  deposited 

sediment.  This  procedure  is  repeated  for  each  daia  set  and  unique  values  of  A„  arc 
found. 

Fukuda  and  Lick  (1980)  report  entrainment  results  for  mud  (67.79t  clay)  from  the 
central  part  of  Lake  Erie  (Table  1).  Sediment  samples  were  allowed  to  consolidate  from 
slurries  lor  periods  of  1  to  10  days.  It  is  assumed  in  this  study  thai  the  sediment  had  noi 
consolidated  significantly  within  one  day,  which  is  equivalent  to  a  water  content  of  74% 
(shaded  rows  in  Table  1).  This  is  the  shortest  consolidation  interval  for  which  data  are 
available.  The  solution  to  Eq.  4  assumes  there  has  been  no  bioiurbation  or  consolidation 
w.c.,  Aa=  Ac=  I )  and  we  are  solving  for  A0  and  m  only.  However,  we  will  use  a  con¬ 
stant  value  of  3  for  the  parameter  m.  following  Lick  ei  «/.  ( 1994).  Fukuda  and  Lick 
(1980)  do  not  report  a  value  for  so  we  use  a  value  of  0.04  Pa.  which  is  reasonable  for 
recently  deposited  mud.  Equation  4  is  fit  to  the  available  high-water-content  data  with 
tinly  Aq  as  an  adjustable  parameter. 


I  able  1 .  Experimental  conditions  and  entrainment  rates  (kg/m/s)  for  Lake  Erie  sedimem 
(Fukuda  and  Lick.  1980).  ^  =  measured;  Em  =  calculated  by  Eq.  4  with  Ac  (initial)  = 
~  L  Aq  =  3.61x10"*  kg/nr/s.  m  -  3,  and  rt  =  0.04  Pa. 


%  (Pa) 

Water 

Content  (%} 

Measured  £$ 

Predicted  Es 

0.0875 

74.0 

8.04x1  O'* 

6.04x10 4 

0.0636 

74.0 

5,94x10* 

7.41xl0’4 

0.0947 

69.5 

1.08xlOJ 

9.23x10 4 

0.0677 

69.5 

1.18x10** 

1. 12x10 4 

0  1091 

64.0 

1 .02x1  O'1 

I.86X10'1 

0.0885 

64,0 

4,79x1  O'6 

6.43x10 4 

0.1133 

61.5 

1.77x10-' 

2.22x10 J 

0.0915 

61.5 

3.70x10-* 

7.70x1  O'4 

0.0616 

61.5 

3.53x1 0*’ 

5.68x10-’ 

£c/E.vi 

(Estimated  At> 

1.3300 
0.8012 
0.1 170 
0,0984 
0.0548 
0.0074 
0.0079 
0.0048 
0.0062 
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The  consolidation  coefficient  A(  for  all  sediment  suites  is  evaluated  using  labora¬ 
tory  entrainment  data  for  the  Lake  Erie  sediment,  which  contained  no  organisms.  The 
value  of  Aq  determined  for  the  high-water-content  data  is  used  to  calculate  entrainment 
rates  for  the  experiments  at  longer  consolidation  intervals,  which  correspond  to  lower 
water  contents  (Table  1 ).  The  ratio  of  the  measured  (Eo)  and  modeled  (£* »)  entrainment 
rates  is  the  estimated  consolidation  coefficient  Ar  =  £o/E*n«  which  is  plotted  against  the 
water  content  in  Pig.  1.  A  function  for  Ac  is  derived  from  these  points.  Note  that  the 
selected  function  is  strongly  dependent  on  the  larger  values  of  Ac  at  high  water  contents. 
The  result  of  this  data-fitting  procedure  will  be  discussed  in  the  results  section. 

An  analogous  procedure  is  used  to  estimate  the  bioturbation  parameter  A8.  We  will 
use  entrainment  data  for  homogeneous  slurries  produced  from  Tamar  Estuary  sediments 
with  infauna  (Lintern  et  ai%  2002).  These  samples  had  consolidation  intervals  of  5  to  58 
days  (Table  2)  and  exhibited  obvious  biological  effects.  Although  these  sediments 
exhibited  an  increase  in  density  with  time,  they  will  be  treated  in  this  study  as  having 
been  dominated  by  bioturbation  (/.£..  Ac  =  1)  as  observed  in  the  entrainment  data,  which 
will  be  discussed  in  the  results  section.  To  determine  a  baseline  coefficient  A0  for  these 
sediments,  we  again  use  samples  with  the  shortest  time  since  deposition,  which  for  these 
sediments  is  5  days  (shaded  rows  in  Table  2),  and  apply  Eq.  4.  Ideally,  the  value  of  An 
can  be  determined  for  a  specific  sediment  type  but,  since  the  mineralogy  of  the  Tamar 


Fig.  1.  Plot  of  the  consolidation  parameter  Ac  for  the  data  from  Table  I.  The  line  is  a 
plot  of  Eq.  5  (r2  =  0.4597).  Note  that  Ac  E  E</E%n  from  Table  1 . 
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I  able  2.  Experimental  conditions  and  entrainment  rates  (kg/m:/s)  for  ISIS  experiments 
(Lintcrn  el  al. .  2002).  £0  =  measured;  EM\  =  calculated  by  Eq.  4  with  Ag  (initial )=  Ac  = 
I  .  4o  =  1.79x10 '  kg/nr/s.  and  m  =  3. 


rh  (Pa) 

T,  (Pa) 

Consolidation 

Interval 

(days) 

Measured  Eq 

Predicted 

£«i 

(Estimated  An) 

0.2 

0.100 

5 

1.97x10'’ 

3.92x10* 

S^xlO1 

0.5 

0.100 

5 

5.44X105 

1.09x1  O’4 

5.0x10'' 

0.1 

0.010 

19 

2.97x10’ 

1.30x10 1 

2.28x10’ 

0.2 

0.010 

19 

4.94x10’ 

1.23x10' 

4.03x10’ 

0.5 

0.010 

19 

1.43x10 4 

2.11 

6.82x10'’ 

0.1 

0.075 

30 

2.11x10’ 

6.63x10 7 

3.18x10' 

0.2 

0.075 

30 

5  19x10’ 

8.29x10’ 

6.27x10' 

0.1 

0.065 

43 

2.97x10’ 

2.79x10* 

1.06x10' 

0.2 

0.065 

43 

9.92x10’ 

1.60x10"* 

6 18x10 ' 

0.1 

0.073 

58 

5.44x10’ 

9.06x1 0’7 

6.01x10' 

0.2 

0.073 

58 

755x10’ 

9.42x10’ 

8.02x10' 

Estuary  mud  is  different  from  that  of  Lake  Eric  mud,  it  is  not  expected  that  Aa  will  be 
tht  same.  I  he  effect  of  mineralogy  is  beyond  the  scope  of  this  work.  The  erosion  data 
from  Lintcrn  et  at.  (2002)  arc  reported  as  total  entrainment  f.  which  is  an  equilibrium 
value  (kg/m*).  This  is  transformed  into  the  entrainment  rate  E0  by  dividing  by  an 
estimated  time  step  of  3600  s.  based  on  the  results  of  Fukuda  and  Lick  (1980).  This 
permits  a  quantitative  comparison  of  the  solutions  of  Eq.  4  for  the  data  contained  in 
Tables  I  and  2. 

The  consolidation  and  bioturbalion  functions  that  are  derived  for  Ac  and  Ag.  re¬ 
spectively.  will  be  evaluated  using  mud  from  Long  Island  Sound  (Tsai  and  Lick.  1988) 
and  undisturbed  cores  from  the  Fox  River  in  Wisconsin  (McNeil  et  al..  1996).  The  Long 
Island  mud  was  processed  to  be  completely  sterile.  Some  samples  were  subsequently 
seeded  with  Nucula,  thus  supplying  an  opportunity  to  examine  the  relative  effects  of 
consolidation  and  bioturbation  on  the  same  sediment.  The  Fox  River  data  is  included  as 
a  demonstration  of  the  many  problems  that  are  encountered  with  cohesive  sediments 
because  of  the  unknown  dcpositional  history,  overall  coarse  gram  size,  and  high  organic 
carbon  content  of  this  sediment 
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4*  RESULTS  AND  DISCUSSION 

This  section  will  use  the  laboratory  entrainment  data  described  above  to  determine 
expressions  for  the  coefficients  in  Eq,  4.  In  all  of  the  work  discussed  below,  the  value  of 
m  is  equal  to  3.  Where  unavailable  from  the  source,  an  erosion  time  step  of  3600  s  will 
be  assumed.  The  value  of  fc  will  be  taken  from  the  sources  where  given*  or  set  to  0.04 
Pa  where  unknown.  We  are  unaware  of  experiments  designed  to  address  this  specific 
problem  so  we  will  use  available  data  with  the  understanding  that  they  may  not  be  well- 
suited  for  our  purposes. 

4*1.  Base  entrainment  coefficient*  At 

The  entrainment  data  with  the  largest  water  content  from  the  Lake  Erie  samples 
(shaded  rows  in  Table  1 )  are  used  to  find  a  best-fit  of  Eq,  4  and  thus  the  value  of  A0 for 
this  suite  of  data.  The  resultant  base  entrainment  of  3;Glxl0^  kg/m2/s  is  used  in  the  next 
section  to  determine  the  appropriate  function  for  Ac -  Similarly,  Eq.  4  is  fit  to  the  Tamar 
mud  with  the  shortest  consolidation  time  and  presumably  least  bioturbaiion  (shaded 
rows  in  Table  2),  and  a  value  of  1.79*10  *  kg/m'/s  Is  estimated  for  A0,  The  functional 
form  of  Ah  will  be  found  below  using  this  value  in  Eq.  4  for  the  data  in  Table  2f  as 
discussed  in  the  previous  section. 

The  Long  Island  Sound  data  are  shown  in  Table  3*  and  the  sediments  with  the  high¬ 
est  water  content  (shaded  rows)  arc  used  to  calculate  for  this  sediment.  The  best  fit  of 
Eq,  4  is  found  for  A a  equal  to  L523xl0'7  kg/irT/s*  with  T*.  =  0  04  Pa  and  a  lime  step  of 
3600  s.  The  sediment  cores  from  the  Fox  River  (Table  4)  were  undisturbed.  A  core  from 
a  water  depth  of  5.9  m  in  Trenton  Channel  consists  of  silt  and  day*  with  gas  bubbles, 
and  surface  oligochucics  and  macrophytes  (McNeil  et  a/..  1996).  This  sediment  has  an 
unknown  history  of  consolidation  and  bioturbaiion  along  with  other  influences  that  are 
not  examined  in  this  study.  In  order  to  determine  A<>  we  need  to  evaluate  the  entrainment 
rate  for  newiy  deposited  sediment*  which  is  unlikely  for  these  samples.  Nevertheless, 
ihis  exercise  is  instructive  in  demoastroling  the  limitations  of  the  model.  We  will  neglect 
bioturbaiion  in  this  sediment  because  we  have  no  information  on  deposition  time, 
whereas  the  water  content  was  measured.  The  following  variables  are  used:  Tf  =  0*1  Pa, 
an  erosion  time  step  of  3600  s,  and  m  =  3,  The  best-fit  value  of  is  2.295x10"  kg/m  Vs. 

The  base  entrainment  rates  for  the  freshwater  sediment  from  Lake  Eric  and  the  Fox 
River  are  much  larger  than  for  the  estuarine  sediment  from  the  Tamar  Estuary  and  Long 
Island  Sound.  This  is  a  result  of  the  salinity  of  seawater,  which  increases  aggregation  of 
cloy  minerals  (Arialhurai  ct  aL ,  1983),  It  should  be  possible  to  identify  a  functional 
form  for  this  process  as  well  ii  appropriate  experimental  data  arc  available.  The  large- 
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labic  V  Experimental  conditions  and  entrainment  rales  (kg/mVs)  for  Long  Island  Sound 
sediment  ( 1  sai  and  Lick*  1988),  £0  is  measured;  £Wi  is  calculated  by  Eq.  4  with  AB  =  Ar 

L  Art  =  i  523x  10  kg/m  /s*  and  m  -  3;  is  computed  the  same  as  £WI  but  with  Ar 
calculated  from  Eq.  5. 


T,.  (Pa) 

Water 

Content 

w 

Measured  £<> 
(kg/rrr/s) 

Predicted  Em 
without 

consolidation 

Predicted 

with 

consolidation 

0  1 

7L3 

5.60x1  O'1 

4,643x!0'7 

1.27x1  O'7  " 

0.2 

71.3 

8.06x10* 

9.74x10* 

2. 67x1  O'7 

0*4 

71 ,3 

1.15x10-* 

UlxIO-* 

3.04x10  s 

0.1 

71 

6  14x10 7 

4.64XI0'7 

1.09x1 0'7 

0.2 

71 

2.61x10* 

1.01x10  s 

2.38x10* 

0.4 

71 

9.5  1x10'* 

1. 08x10"* 

2.55x10* 

0,8 

71 

1.41x10-* 

1.04x10-* 

2.47x10"' 

0*2 

69.9 

1.34x10* 

9.74x1 0‘* 

1.32x10* 

0.4 

69.9 

2. 90x1  O'* 

1.11x10* 

1.50x10* 

0.8 

69.9 

9,83x10'* 

1.04x10'* 

1.41x10* 

1.2 

69.9 

3.73x10** 

9.56x10'* 

1.29x10-* 

I  able  4,  Entrainment  rate  for  Fox  River  sediments  (McNeil  et  at..  1996),  The  model 
predictions  <£Af  1 )  were  calculated  using  A<>  =  2,295* I  O'1  kg/m2/s*  Ag=  L  Ac  =  Eq.  5,  ic 
^OJPa,  an  equilibrium  time  of  3600  s,  and  m  -  3. 


Depth 

(m) 

n  (Pa) 

Water  Content 

<%) 

Measured  Eg 
(kg/nr/s) 

Predicted 

(kg/nr/s) 

0,017 

0,25 

64*8 

7.74 

8.13x10-* 

0,27 

0.25 

50.0 

4x10-* 

3.07x10* 

0*26 

0.6 

50,0 

5.28x10 J 

1.14X10'5 

0.34 

LI 

60.0 

6.36xl02 

9.1  lxl0': 

0,15 

2.2 

49.0 

8.57x10  J 

8.44X10'1 

0.22 

4.5 

49,0 

4.18xl0l 

7.76 
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base  entrainment  rale  for  l he  Fox  River  data  further  suggests  that  bioiurhation  may  also 
be  an  important  process  at  this  location,  as  indicated  by  the  presence  of  organisms  at  the 
river  bed, 

4.2.  The  consolidation  coefficient,  Ac 

The  entrainment  data  for  the  Lake  Erie  sedtmenl  are  used  to  determine  an  expres¬ 
sion  for  Ar  because  they  have  the  largest  range  of  water  content  (6L5  -  74%),  and 
presumably  consolidation.  The  values  of  Ed  Em  from  Table  1  are  an  estimate  of  the 
consolidation  coefficient  A<%  which  is  plotted  against  water  content  in  Fig  I.  An  expo¬ 
nential  function  is  fit  to  these  data: 

AfsC|e|H''*#l"  (5) 

where:  C,  -  0.00397.  I  =  1.99315*  W  -  water  content  (%},  and  WQ  =  62.85  for  this 
study.  The  squared  correlation  coefficient  r2  for  Eq.  5.  as  plotted  in  Fig.  i.  is  0,4597. 
The  range  of  water  content  for  which  Eq.  5  is  valid  is  limited  by  the  data  used  to  derive 
it.  which  in  this  case  is  61.5  to  74%.  An  exponential  function  is  reasonable  because 
consolidation  can  be  represented  to  first-order  as  an  exponential  process.  However,  the 
constants  Cj,  and  t  arc  dependent  on  the  specific  data  used  to  find  Ao  and  their 
values  will  change  as  other  data  are  examined  as  will  the  functional  form  of  Eq,  5. 
Equation  5  implies  that  as  sediment  is  compacted  and  W  decreases,  the  entrainment  rate 
will  decrease  as  well.  With  respect  to  the  model,  this  indicates  that  small  values  of  Ac 
correct  for  the  base  entrainment  rate  over- predicting  the  entrainment,  because  of  an 
increase  in  Tf.  This  error  increases  with  consolidation.  If  the  data  used  to  calculate 
are  properly  selected,  Ac  should  always  be  less  than  unity  but  this  is  unlikely  due  to 
measurement  uncertainties  and  the  influence  of  processes  not  included  in  the  model. 

4.3.  The  hint ur billion  coefficient,  A # 

The  laboratory  data  for  slurries  from  the  Tamar  Estuary  are  used  to  examine  the 
functional  form  of  the  bioturbation  coefficient  A *  because  they  are  dominated  by  bio¬ 
logical  processes  and  the  results  are  well  documented.  Before  determining  the  biolurba- 
lion  function,  it  is  useful  to  note  the  difficulty  in  using  Eq.  4  for  biologically  influenced 
sediment.  The  ratio  Ed  Em  for  the  low-  shear- stress  data  used  to  calculate  (top  row  in 
Table  2)  exceeds  500.  which  is  significantly  greater  than  any  other  value  in  the  table. 
This  occurs  because  the  (i/^  -  1)  term  in  Eq,  4  is  less  than  unity  for  this  experiment, 
resulting  in  a  substantial  decrease  in  the  predicted  entrainment.  This  may  be  due  to 
many  causes  that  are  not  included  in  the  model.  This  phenomenon  has  been  noted  by 
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Lavcllc  et  ai  (1984)  and  it  is  of  continuing  interest  in  predicting  cohesive  sediment 
entrainment.  We  will  not  use  the  error  data  from  this  experiment  in  deriving  AB%  how¬ 
ever.  because  it  skews  the  resulting  function. 

1  he  bioturbation  parameter  AB%  which  is  estimated  from  the  ratio  of  £<>/£w ,  (Table 
2).  is  plotted  as  a  function  of  consolidation  time  in  Fig.  2.  These  data  are  represented  by 
the  following  expression: 

Ah  =C2  +  Cytd  (6) 

where:  C2  =  5.093x10  s;  C,  =  0.02186;  C4  =  6.85x10-*;  and  /,  =  consolidation  time 
(days).  The  squared  correlation  coefficient  r2  is  0.3563,  which  is  rather  low  because  of 
the  small  sample  size  and  range  of  reported  entrainment  rates.  The  constants  C2,  C3.  and 
C4  will  have  different  values  if  other  experimental  data  are  used  to  find  AB.  but  we 
propose  that  Eq.  6  is  a  reasonable  function.  These  data  do  not  fit  an  exponential 
function  as  well  as  the  consolidation  data  because  of  the  complex  impacts  of  biological 
processes.  For  example.  I-intcm  et  al.  (2002)  noted  the  presence  of  diatoms  and  other 
biological  material  w'ithin  the  homogenized  sediment,  in  addition  to  complex  structures 
caused  by  burrowing  organisms.  Furthermore,  not  all  biological  activity  increases 
credibility;  for  example,  biofilms  promote  stability  in  sediments  whereas  bioturbation 
destabilizes  the  bed.  Overall,  these  deposits  were  mostly  bioturbated.  which  is  the 


Bioturbation  Time  (Days) 

Fig  2.  Plot  of  the  bioturbation  coefficient  Ati  for  the  data  from  Table  2.  The  line  is  a 
plot  of  Fq.  6  (  r  =  0.3563).  The  coefficient  AB  =  £</£%//  from  Table  2.  Note  that  the 
first  data  entry  is  not  plotted  See  text  for  explanation. 
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reason  Fig.  2  shows  a  slight  increase  in  entrainment  with  time.  If  the  data  used  to 
calculate  Ao  are  representative  of  newly  deposited  sediment,  AH  should  always  be 
greater  than  unity  because  of  the  increasing  disruption  as  organisms  hioturbatc  ihe  bed. 
However,  as  seen  in  Fig.  2f  A&  is  less  than  l  for  most  bioturbation  times  because  of  the 
experimental  uncertainly  (note  the  3  data  points  at  19  days). 

4.4,  Comparison  to  other  entrainment  data 

The  sediment  samples  studied  by  Tsai  and  Lick  (1988)  were  collected  in  a  water 
depth  of  20  m  in  Long  Island  Sound.  They  contained  dark  grey  organic  sill  and  clay 
with  959f  finer  than  74  pm.  The  samples  were  homogenized  and  placed  in  an  annular 
flume  0.204  m  deep.  The  water  content  at  the  time  of  testing  varied  from  70.4  to  72.2*fr. 
This  is  a  small  range  and,  in  fact,  there  was  some  overlap  in  water  content  during  the 
consolidation  interval.  The  value  of  te  was  not  reported  so  we  use  0.04  Pa.  A  compari¬ 
son  will  be  made  with  the  first-hour  resuspension  rate  (Fig.  4  from  Tsai  and  Lick.  1988), 
so  a  time  step  of  3600  s  is  used.  The  value  of  Ao  (1.523x10  kg/mVs)  was  determined 
from  samples  that  were  measured  after  I  day  of  consolidation  as  discussed  above. 

The  entrainment  data  for  these  samples  arc  given  in  Table  3.  There  is  a  significant 
trend  of  decreasing  entrainment  rate  with  lime  (correlated  with  water  content)  for  the 
observations  at  all  shear  stresses  (open  squares  in  Fig.  3),  This  behavior  is  what  wc 
would  expect  for  consolidation.  Note  that  the  consolidation  interval  labels  l  D)  on 
Fig.  3  refer  to  the  measurements  only.  Equation  4  is  solved  using  the  base  entrainment 
iAfl  =  At  -  I )  and  plotted  on  Fig,  3  as  well  (open  circles.  Em  in  Table  3).  Note  that  the 
data  points  for  all  water  contents  fall  upon  each  other,  because  of  the  lack  of  consol ida- 
lion  when  Ac  =  l  This  is  the  reason  it  is  necessary  to  adjust  other  parameters  for  Eqs.  I 

-  3.  The  prediction  is  not  bad  because  the  total  consolidation  interval  was  only  7  days. 
The  effects  of  consolidation  are  seen  at  ail  shear  stresses  when  Ac  is  computed  using  Eq. 
3,  however.  The  results  (crosses  in  Fig.  3)  match  the  measurements  better  as  well  (Ew: 
in  Table  3). 

Some  of  the  homogenized  sediment  from  Long  Island  Sound  was  seeded  with  de¬ 
posit-feeding  Nucuta  clams  at  a  density  of  1 15CVuf  (Tsai  and  Lick,  1988).  This  signifi* 
cantlv  increased  the  equilibrium  concentration  of  sediment  at  all  shear  stresses,  Tsai  and 
Lick  ( 1988)  report  the  first-hour  entrainment  rate  for  the  sterile  data  but  only  the  mean 
concentration  C*f  at  different  shear  stresses  for  the  seeded  samples.  Nevertheless,  the 
entrainment  rate  can  be  approximated  from  CM  by:  E  -  -  Cju(i)1/M  where  h 

-  the  flume  depth  (0.204  m)  and  At  =  the  time  to  reach  CM  for  an  increase  in  ih.  This 
calculation  was  performed  on  the  data  of  Tsai  and  Lick  t their  Fig.  8)  for  comparison 
with  the  values  predicted  by  Eq.  4  using  Eq.  6  to  find  the  bioturbation  coelficient  A^ 
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l  ig  3  Plot  of  the  entrainment  data  for  Long  Island  Sound  sediment:  measured  (Eq) 
(Tsai  and  Lick,  1988)  (squares);  calculated  by  Eq.  4  with  AB  ^  Ac  =  I  {Em,  pluses); 
and  Eq.  4  with  Ac  calculated  by  Eq.  5  and  Ac  =  1  ( E ^  Circles)*  The  consolidation 
intervals  in  days  for  the  measured  data  are  indicated. 


1  he  results  (Fig.  4)  show  that  Eq.  6  is  not  general  enough  for  wide  application  because 
it  was  derived  from  a  limited  database  with  somewhat  unique  characteristics* 

The  entrainment  rates  for  the  undisturbed  sediment  from  the  Fox  River  were  meas¬ 
ured  using  SEDFLUME  (McNeil  ei  a/*,  1996).  This  sediment  contained  silt  and  clay 
(mean  grain  size  less  than  20  pm),  gas  bubbles,  infauna,  grass,  and  organic  carbon 
(approximately  10%  i.  The  entrainment  rate  was  measured  at  a  range  of  depths  within 
the  core  for  shear  stresses  ranging  from  0,25  -  4.5  Pa.  The  water  content  varied  from 
64.8%  at  the  surface  to  48.1%  at  depth,  indicating  that  this  sediment  was  not  freshly 
deposited,  A  base  entrainment  rate  (4q)  of  2,259x10  '  kg/rrr/s,  is  calculated  as  described 
m  section  4.1  and  Ac  is  computed  using  Eq.  5.  The  bulk  density  required  to  calculate  E 
is  found  using  bulk  density  data  from  McNeil  ei  aL  (I9%#  Fig.  8);  the  surface  bulk 
density  of  the  core  was  1 166  kg/m\  The  entrainment  rates  for  the  measured  samples 
(squares  in  Fig.  5)  vary  at  smaller  values  of  and  increase  with  shear  stress.  The 
entrainment  rates  predicted  by  Eq,  4  increase  with  shear  stress  also  while  being  an 
order-o t- magnitude  high  at  large  .shear  stresses.  The  predicted  entrainment  rates  are 
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Shear  Stress  (Pa) 


Fig.  4.  Plot  of  the  entrainment  data  for  Long  Island  Sound  sediment  seeded  with 
Nucula  after  7  days  (Tsai  and  Lick,  1988).  The  following  conditions  were  used  for  the 
simulation:  time  step  =  3600  s;  ic  =  0  04  Pa;  A calculated  from  Eq.  6;  Ac  =  I;  Aq  = 
1.523x10 7  kg/m2/s. 


probably  too  high  because  of  post-dcpositional  modification  of  the  sediment.  The 
consolidation  expression  (Eq.  5)  is  based  on  limited  data  from  deposited  slurries. 
Despite  the  difficulties  of  predicting  the  behavior  of  this  sediment,  the  fit  is  reasonable, 
especially  since  no  curve-fitting  process  was  used  for  consolidation.  It  should  also  be 
noted  that  the  base  entrainment  coefficient  for  this  sediment,  which  was  estimated  from 
a  surface  sample,  probably  does  not  incorporate  the  effects  of  burial  as  deep  as  0.4  m. 


5.  CONCLUSIONS 

Consolidation  and  bioturbation  appear  to  have  opposing  effects  on  cohesive  sedi¬ 
ment  entrainment  for  most  of  the  sediment  studied  herein.  For  example,  bioturbation  by 
infauna  reworks  sediment  and  increases  erodibility  for  these  samples  whereas  consolida¬ 
tion  and  dewatering  tends  to  decrease  entrainment  rates.  Ihus.  these  preliminary  results 
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I  ig.  5.  Entrainment  rates  for  Fox  River  sediments  (McNeil  et  aL,  1996).  The  following 
conditions  were  used  for  the  simulation:  time  step  =  3600  s;  %  =  0.1  Pa;  A,-  calculated 
from  Eq.  5M#  =  I ;  ,\0  =  2.259x10  *  kg/nr/s. 


suggest  that,  as  a  first  approximation,  these  processes  can  be  treated  as  mutually  exclu¬ 
sive.  It  bioturhation  is  significant,  consolidation  can  probably  be  neglected  and  Af 
assumed  to  be  t.  whereas  if  bioturbation  is  negligible  (i.e.,  rapid  deposition  and/or  a 
hypoxic  water  column),  the  sediment's  physical  properties  are  probably  dominated  by 
consolidation  or  other  processes  that  decrease  credibility  (e.g..  organic  content).  It  is 
not  suggested  that  this  simple  approach  incorporates  all  environmental  effects  but  it  is 
usable  until  better  entrainment  data  become  available.  The  greatest  advantage  of  the 
modular  entrainment  model  (Eq.  4)  is  that  it  permits  us  to  examine  this  complex  process 

with  a  more  physically  realistic  approach  than  parameter  calibration  based  on  entrain- 
mcm  daia. 
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